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Abstract. This article gives an overview of many of the recent developments in 
understanding the structure of relativistic scattering amplitudes in gauge theories 
ranging from QCD to Af = 4 super- Yang-Mills theory, as well as (super)gravity. I 
also provide a pedagogical introduction to some of the basic tools used to organize 
and illuminate the color and kinematic structure of amplitudes. This article is an 
invited review introducing a special issue of Journal of Physics A devoted to "Scattering 
Amplitudes in Gauge Theories" . 
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1. Overview 

The astrophysicist Subrahmanyan Chandrasekhar once referred to black holes as "the 
most perfect macroscopic objects there are in the universe: the only elements in their 
construction are our concepts of space and time." p] This special issue of Journal of 
Physics A is devoted to a completely different topic, the quantum-mechanical amplitudes 
for the scattering of relativistic particles, in which a revolution in our understanding has 
been taking place recently However, scattering amplitudes represent a kind of short- 
distance complement to the perfection Chandrasekhar saw in black holes. We might 
consider them to be the most perfect microscopic structures in the universe. 

Black holes are gigantic and (almost) eternal. Scattering amplitudes, in contrast, 
portray events that wink in and out of existence much faster than the blink of an 
eye. They describe processes at the shortest distance scales that can be probed in the 
laboratory. When gravitons scatter, interacting via Einstein's field equations, again in 
some sense the only elements in the construction are our concepts of space and time - 
although the boundary conditions at infinity can be more detailed than for black holes, 
involving multiple plane-wave ripples (the gravitons) propagating in various directions 
through Minkowski space-time. When other relativistic particles scatter, as in non- 
abelian gauge theory, more information has to be provided, such as the gauge group, the 
coupling constant, the matter content of the theory, the particle masses, and the precise 
set of particle types and spins being scattered. However, in the ultra-relativistic limit, 
if there are only gauge interactions, much of the group-theoretical information can be 
stripped from the amplitudes, and certain color-ordered "primitive" amplitudes can be 
defined which are quite universal. These amplitudes have intricate analytic properties, 
which have been teased out over the decades, but especially in the past few years. A 
variety of new, efficient techniques have been found to construct amplitudes, and hidden 
symmetries have been identified which illustrate the amplitudes' "perfection." 

In the real universe, black holes are not found in an absolute vacuum. Typically 
they are surrounded by accretion disks of gas and dust, which muddy their perfection 
somewhat, but also enrich the range of phenomena involving them, and make them 
much easier to observe, albeit indirectly. Similarly, scattering amplitudes for quantum 
chromodynamics, or QCD, the most relativistic non-abelian gauge theory we can study 
in the laboratory, are shrouded from direct view by confinement. The quarks and gluons 
that theorists scatter mentally are never seen experimentally. In the initial state of a 
hadronic collision they are bound into a proton or other hadron. In the final state they 
emerge as collimated jets of particles. Despite this fact, quantities that are sensitive 
only to very short distances, so-called infrared-safe observables such as jet production 
cross sections, can be computed reliably in terms of quarks and gluons, in a systematic 
expansion in the strong coupling a s . 

If we could go to asymptotically high energies, so that ct s were infinitesimal, it would 
be enough to consider just the leading order in this expansion, which is generated by the 
Born approximation, or tree-level amplitudes. The tree-level amplitudes of QCD have 
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even more "perfection" than the generic, loop-level amplitudes, because they coincide 
with the tree amplitudes of a much more symmetric theory, M = 4 supersymmetric 
Yang-Mills theory (TV = 4 sYM). This theory, which is the subject of many of the 
articles in this special issue, has the maximum amount of supersymmetry possible in 
a gauge theory. Its four supercharges can be used to transform the helicity +1 gluon 
state, by \ unit of helicity at a time, all the way to the helicity —1 CPT conjugate 
gluon state. All states in the theory belong to a single supermultiplet, transforming 
in the adjoint representation of the gauge group. Therefore all interactions are related 
by supersymmetry to the triple-gluon vertex, and the theory has a single dimensionless 
coupling g. The j3 function for M = 4 sYM vanishes exactly for any gauge group, 
(3(g) = 0, so that the theory is exactly conformally invariant at the quantum level, as 
well as classically. (The classical conformal invariance of QCD with massless quarks is 
of course spoiled by its nonvanishing one-loop (3 function, which leads to asymptotic 
freedom.) 

In investigations of scattering amplitudes in M = 4 sYM, the gauge group under 
consideration is usually SU(iV c ). Quite often, the limit N c — > oo is also taken, at a fixed 
value of the 't Hooft coupling, A = g 2 N c . For A 1, perturbation theory can be applied, 
and for large N c only planar Feynman diagrams contribute. For A 3> 1, perturbation 
theory breaks down, but the AdS/CFT correspondence [2] can be used. In the large- 
N c or planar limit of M = 4 sYM, even more remarkable symmetries and relations 
emerge for scattering amplitudes. As we will discuss further below, these properties 
include dual conformal (super) symmetry, which is part of a larger Yangian invariance; 
a description of amplitudes in terms of spaces of complex planes (Grassmannians); and 
in terms of various types of twistors; and a relation between scattering amplitudes and 
the expectation values of Wilson loops for closed polygons bounded by light-like edges. 
Many of these properties are intimately connected with a strong-coupling picture of 
gluon scattering in terms of strings moving in anti-de-Sitter space [3]. 

The exceptional simplicity and numerous hidden symmetries of M = 4 sYM have 
made it a playground for theorists interested in scattering amplitudes for a couple of 
decades. This interest has accelerated rapidly in the past few years. Because the tree- 
level amplitudes of QCD and A^ = 4 sYM coincide, the early discoveries about QCD 
helicity amplitudes, motivated by the physics of jet production, were also discoveries 
about M = 4 sYM. For example, the Parke- Taylor formula jl] for the maximal-helicity- 
violating (MHV) sequence of n-gluon amplitudes, which is just a single-term expression, 
even for an arbitrarily large value of n, is also consistent with all the symmetries ofAf = 4 
sYM [5], including some that were only unveiled decades later. 

Many general properties of tree amplitudes were understood early on, such as their 
factorization properties [6] and the (A/" = 1 or M = 2) supersymmetry Ward identities 
that they obey [TJ E\- However, a fuller understanding and exploitation of these 
properties have come only in recent years. For instance, the complete solutions of 
the long-known Af = 4 supersymmetry Ward identities, and the corresponding M = 8 
identities in supergravity, have been worked out recently by Elvang, Freedman and 



Introduction to scattering amplitudes 



4 



Kiermaier [5J, and are reviewed here |1U] . 

Several years ago, Witten [11] discovered that the twistor space developed decades 
earlier by Penrose [12] gave a natural description of tree-level amplitudes in massless 
gauge theory. As discussed in section [3j such amplitudes have a natural description in 
terms of two-component Weyl spinors, a left-handed and a right-handed spinor for each 
external state. To go to twistor space, one performs a Fourier transform on half the 
variables, namely the left-handed spinors, leaving the right-handed ones alone. Witten 
found that gauge theory could be reformulated in terms of a topological string moving in 
twistor space. Gauge amplitudes turn out to be localized on particular curves in twistor 
space. In one approach, the curves are intersecting lines: a single line for the simplest 
MHV amplitudes, a pair of lines for the next-to-maximally-helicity-violating (NMHV) 
amplitudes, and so on. This version gave rise to the CSW or MHV rules for gauge 
theory |13j . whose MHV vertices are a particular off-shell continuation of the Parke- 
Taylor formulae. These developments are reviewed here by Brandhuber, Spence and 
Travaglini [H], and by Adamo, Bullimore, Mason and Skinner [15]. The latter reference 
also discusses another class of twistors, namely momentum twistors, which also provides 
a natural set of kinematic variables for amplitudes, in that they automatically satisfy the 
kinematic constraints of momentum conservation and the mass-shell conditions. Finally, 
the correspondence between Wilson loops and scattering amplitudes is described in this 
article [15] from the perspective of twistor space. 

Britto, Cachazo, Feng and Witten (BCFW) [16] recognized that factorization 
properties are powerful enough to allow the derivation of recursion relations [IT] for 
tree amplitudes, in terms of on-shell lower-point amplitudes, which are evaluated at 
particular complex momenta. This observation, reviewed here by Brandhuber, Spence 
and Travaglini [14], makes essential use of the analyticity, or plasticity, of amplitudes; 
that is, how they vary under smooth deformations of the kinematics. BCFW embedded 
an amplitude A = A(0) into a one-complex-parameter family of on-shell amplitudes 
A(z) with shifted complex momenta, and associated the poles of A(z) in the z plane 
with factorization of the amplitude onto simpler lower-point amplitudes. One can 
also construct a supersymmetrized version of the BCFW recursion relation, in which 
Grassmann variables r] A associated with J\f = 4 supersymmetry are shifted along with 
the momenta [T8l fl~9| IT4"] . This recursion relation can be solved explicitly to yield all 
tree amplitudes in M = 4 sYM [20]; the solution and its many symmetries are reviewed 
here by Drummond [2T] . 

Scattering amplitudes in massless gauge theories all have quite similar structure 
at tree level, once the color factors have been stripped off, as discussed in section [2j 
On the other hand, the structure of loop amplitudes depends critically on the theory. 
The simplicity of one-loop multi-gluon amplitudes in M = 4 sYM played a key role 
in the development of the unitarity method [22J, [23]. This method reconstructs loop 
amplitudes from their unitarity cuts, again exploiting analyticity. In general at one 
loop, as explained in this issue by Britto [21] and Ita [25], unitarity cuts are matched 
against a decomposition of the amplitude in terms of a set of scalar integrals — boxes, 
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triangles, bubbles and (sometimes) tadpoles — in order to determine the coefficients 
of the integrals. In M = 4 sYM, the high degree of supersymmetry implies that only 
boxes have non-vanishing coefficients. Using unitarity, an infinite sequence of one-loop 
amplitudes (the MHV amplitudes) could be determined in M = 4 sYM from just the 
product of two tree-level MHV (Parke- Taylor) amplitudes [22] . 

Later, generalized unitarity was applied to the problem of computing one-loop 
amplitudes, by extracting additional information from products of three and four tree 
amplitudes, the so-called triple (26] and quadruple (27] cuts. Quadruple cuts have 
four on-shell conditions. These four equations determine the four components of the 
(four-dimensional) loop momenta completely, up to a (possible) two-fold degeneracy. 
This realization allows the box coefficients to be computed very simply in terms of 
the product of four tree amplitudes, glued together as the four corners of a box, and 
evaluated at each of the two solutions for the loop momenta [27]. For M = 4 sYM, there 
are no other coefficients to determine, so the one-loop problem has been "reduced to 
quadrature." (It turns out that maximal M = 8 supergravity can be proven to have the 
same "no-triangle" property [281 EE] , so it too is solved at one loop by the quadruple 
cuts.) 

Generalized unitarity can also be applied very effectively at the multi-loop level. In 
principle it can be used for any gauge (or gravitational) theory. As an example, the two- 
loop four-gluon scattering amplitudes in QCD have been computed in this way [29]. In 
practice the method has been pushed the furthest in M = 4 sYM. The basic techniques 
of multi-loop generalized unitarity are reviewed here by Bern and Huang [30], and by 
Carrasco and Johansson |31j . 

Returning now to one-loop amplitudes for generic gauge theories, including QCD, 
the triangle and bubble coefficients have to be determined, as well as certain rational 
parts, which have no unitarity cuts in four dimensions. The triple cuts contain 
information about the triangle coefficients, but they also receive contributions from 
boxes. Similarly, the ordinary double cuts determine the bubble coefficients, once 
the contributions of triangles and boxes are removed. This separation can be done 
analytically, by making use of the different analytic behavior of the different types of 
cut integrals, as reviewed by Britto [24]. It can also be done numerically, by a suitable 
sampling of the continuum of loop momenta solving the triple or double cut on-shell 
conditions [321 [331 [M] , as described by Ita [25]. These two articles also review methods 
for computing the rational parts of one-loop amplitudes. Although the rational parts 
have no unitarity cuts in four dimensions (by definition), they can be constructed via 
their cuts in dimensional regularization in D = 4 — 2e dimensions. Alternatively, it is 
possible to construct them recursively in the number of legs [34] . 

On-shell loop amplitudes in massless non-abelian gauge theories always contain 
infrared divergences, due to the exchange of soft gluons or virtual collinear splittings. In 
QCD it is conventional to regulate these divergences, as well as the ultraviolet ones, using 
dimensional regularization. In supersymmetric theories, dimensional reduction |35j, 
or the related four-dimensional helicity scheme |36j . can be used to keep the number 
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Figure 1. Factorization of soft and collinear singularities. 

of bosonic and fermionic states the same, and preserve Ward identities for ordinary 
supersymmetry. 

The general structure of the infrared divergences is well understood from decades 
of work in QED as well as QCD [37]. It has been worked out in the most detail in the 
context of dimensional regularization [38] [39J HP] . The basic picture is shown in fig. [TJ 
Soft divergences and collinear divergences associated with the amplitude M each have a 
universal form. They can be factorized from each other and from a hard, short- distance 
part of the amplitude. Soft divergences, denoted by the blob marked S, come from 
exchange of long- wavelength gluons. These gluons do not have the resolving power to 
probe the internal structure of the "jets" J of virtual collinear particles which capture 
the collinear divergences. Soft gluons can only see the overall color charges of the jets. 
An individual jet function Jj depends on the type of particle i, but not on the full 
amplitude kinematics. The soft function S does not depend on the particle types, but 
only on their momenta and color quantum numbers. In general these quantum numbers 
can be mixed by gluon exchange, so S is a matrix in color space. The hard function h 
has no infrared singularities, but generically depends on the particle types, colors and 
kinematics. 

In the planar or large- N c limit, the picture simplifies considerably, to that shown 
in fig. [2l Now M represents the coefficient of a particular color structure, such as 
Tr (T ai T a2 ■ ■ -T an ) (assuming that all external states are in the adjoint representation; 
see section [2]). In the planar limit, individual soft gluons can only connect color- adjacent 
external partons. There is no mixing of different color structures at large N c . One can 
absorb the entire soft function S into jet functions, which corresponds to breaking up 
the right-hand side of fig. [2] into n wedges. Each wedge is bounded by two hard lines, 
and is composed of "half" of each of the two jet functions, as well as the soft gluons 
exchanged between them. Up to nonsingular terms, the wedge controlling the infrared 
divergences represents the square root of the Sudakov form factor, which is defined as 
the amplitude for a color-singlet state to decay into a pair of (adjoint) gluons. 

In dimensional regularization, the Sudakov form factor obeys a particular differ- 
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Figure 2. Soft-collinear factorization in the planar limit. 

ential equation |38j, whose source term is called the cusp anomalous dimension [41]. 
There is an additional constant of integration, called Qq. These two functions (along 
with the function in a non-conformal theory) control the infrared divergences to all 
orders for any amplitude in a planar massless gauge theory. In planar J\f = 4 sYM , 
the cusp anomalous dimension has been determined to all orders in the coupling using 
integrability [32] ■ The BDS ansatz [32] was built up from this description of the infrared 
singularities of planar amplitudes, plus the observation for the two- and three-loop four- 
point amplitudes that the hard function h essentially reduces to a constant, independent 
of the kinematics. 

Besides the cusp anomalous dimension, another place that integrability certainly 
enters multi-loop scattering amplitudes is in multi-Regge-kinematics, a particular class 
of high-energy or small-angle limits of 2 — > (n — 2) particle scatterings. These kinematics 
gave the first indication that the BDS ansatz for MHV amplitudes has to be corrected 
beginning at the six-point level [33] • (It was previously argued [32] using the properties 
of Wilson loops, that the ansatz should be corrected at two loops for a large value 
of n, unless the correspondence between MHV amplitudes and Wilson loops were to 
break down. The hexagon Wilson loop was then computed and found to differ from 
the ansatz prediction [36].) The dynamics of gluons in the transverse plane is given, in 
the leading-logarithmic approximation, by the Hamiltonian for a integrable open spin 
chain, as reviewed here by Bartels, Lipatov and Prygarin |47j . 

In planar TV = 4 sYM, another infrared regulator is more convenient for many 
purposes, in particular for exploring the loop-level consequences of dual conformal 
symmetry, which is not preserved by dimensional regularization. The features of a 
recently-developed "Higgs regulator" [48J are reviewed in this issue by Henn [39]. For 
this regulator, vacuum expectation values are given to some of the adjoint scalar fields 
in the theory, breaking the gauge symmetry in such a way that (in the planar limit) 
the amplitudes are regulated by massive particles circulating around the outside of the 
loop diagrams. Dual conformal symmetry remains exact in a certain sense. When 
the particle masses m, are taken to be much less than the momentum-invariants for the 
scattering process, logarithmic divergences develop, the analogs of the 1/e infrared poles 
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Figure 3. Hierarchy of simplicity in scattering amplitudes for various types of 
gauge theory. 



in dimensional regularization. 

There is a hierarchy of simplicity in the scattering amplitudes for various types of 
gauge theory, as sketched in fig. El This hierarchy begins to be revealed at one loop. The 
outer region of the diagram stands for a generic gauge theory with massive matter fields, 
and perhaps massive gauge bosons, if the gauge symmetry is spontaneously broken, as 
in electroweak theory. One-loop amplitudes in such a theory generically contain tadpole 
integrals. One-particle cuts are nontrivial, and are particularly delicate because of 
external-leg contributions [50j[2l]. The cut structure of loop integrals containing massive 
propagators in the loop is generically somewhat more complicated than the purely 
massless case. Massive particles in the loop can be unstable, which usually necessitates 
complex masses. When one enters the "massless" ring in fig. [3j corresponding to massless 
gauge bosons and matter fields, most of these complications vanish, although there are 
still generically rational parts to compute. The ring "sYM" stands for supersymmetric 
gauge theories. Their one-loop amplitudes can be constructed from four-dimensional 
unitarity cuts alone, i.e. there are no non-trivial rational parts |23j. 

Moving further inward in fig. [3j we arrive at M = 4 sYM. As mentioned earlier, 
at one loop the coefficients of bubble and triangle integrals now vanish, as well as the 
independent rational parts. (There are other gauge theories with vanishing bubble and 
triangle coefficients, at least for their n-gluon amplitudes [ED, El].) The theory becomes 
conformally invariant. It has been conjectured that the leading singularities — the 
multi-loop analogs of the quadruple cuts — are sufficient to determine the amplitudes 
at any loop order [18]. In addition, scattering amplitudes have empirically a predictable, 
uniform transcendental weight [52l [53] . This weight refers to their construction out of 
polylogarithms, logarithms, and Riemann ((n) values. For example, the finite (O(e )) 
terms in one-loop Af = 4 sYM amplitudes are of weight two: They contain some terms 
proportional to the polylogarithm Li2, and others which are products of two logarithms, 
or proportional to C(2), but they do not contain any terms of lower transcendentality. 
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(At one loop this result just follows from the absence of bubble and tadpole integrals.) 
Wherever analytic results are available, this uniform transcendentality property holds. 
For example, at two loops (weight four for the finite terms) it holds for the non-planar, 
subleading-color terms as well as the planar ones |54j. 

Finally, for planar M = 4 sYM, as mentioned earlier, additional symmetries appear, 
dual superconformal invariance and an associated Yangian symmetry. The extra 
symmetries are related to integrability. They give rise to the prospect of adapting 
methods developed for other integrable systems to determine the S"-matrix exactly, 
at least in the large N c limit. Indeed, as reviewed by Drummond [21], an integrand 
exhibiting manifest Yangian invariance has been proposed recently [55] for the general 
L-loop n-point amplitude, based on the BCFW recursion relations; see also ref. [56J. 




At the level of integrated amplitudes, dual conformal invariance has an anomaly 
due to the need for an infrared regulator. The anomaly was first understood in terms 
of Wilson loops rather than amplitudes [57] (for which divergences are ultraviolet in 
nature, rather than infrared). Nevertheless, dual conformal symmetry completely fixes 
the form of the four- and five-point amplitudes. The form to which they are fixed turns 
out to be precisely the ABDK/BDS ansatz [581 S3], which was developed earlier based 
on the structure of infrared divergences of amplitudes, plus patterns observed at two 
and three loops. As discussed in section I3.3[ for six external legs one can build three 
different combinations of kinematic variables that are invariant under all dual conformal 
transformations, the cross-ratios Ui, «2 and u%. At the six-point level, the ABDK/BDS 
ansatz fails at the first nontrivial order, two loops [59] . and dual conformal symmetry 
allows for an additional "remainder function", Rq(ui, 112, %), which depends on the 
three cross-ratios. 

At strong coupling, large A and large- N c , Alday and Maldacena used the AdS/CFT 
correspondence to map the problem of gluon scattering to that of strings moving in anti- 
de-Sitter space [3]. The AdS background is weakly curved in this limit, relative to the 
string scale. Another way of saying this is that the two-dimensional string world-sheet 
stretches over a large area, relative to the scale on which strings fluctuate, allowing a 
semi-classical expansion to be used at large A. In the leading term in the expansion, 
different helicity configurations are not distinguished. The amplitude is given simply by 
exp(— S c \), where S c \ is the classical action that minimizes the area (at least for scattering 
configurations with a Euclidean interpretation). Using a T-duality transformation 
in string theory, Alday and Maldacena showed that the boundary conditions for the 
world-sheet are closed polygons, where each edge is a light-like segment given by the 
momentum ki of the 2 th gluon. The solution to the minimal area problem was found 
explicitly for four-gluon scattering. More recently, in order to handle more complicated 
kinematical configurations than the four-point case, the integrability of the string sigma- 
model action [61)] has been exploited. The minimal area problem has been solved by 
mapping it to a system of equations identical to those of the Thermodynamical Bethe 
Ansatz |61j . 

The first computation at strong coupling was for the four-gluon amplitude, 
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and matched precisely the prediction of the ABDK/BDS ansatz [3]. But it also 
suggested a weak-coupling correspondence between scattering amplitudes and Wilson 
loops for closed light-like polygons, which was rapidly established for a variety of 
cases, initially for the simplest case of MHV scattering amplitudes [62j [57J EH]- The 
Wilson-loop correspondence has been very important both conceptually and technically. 
Conceptually, it is still not entirely clear why it happens at all at weak coupling. 
Technically, given its existence, it allows amplitudes to be computed in terms of 
Wilson line integrals, which are generally somewhat more tractable than Feynman 
loop integrals. Some recent conceptual advances have come from considering certain 
correlation functions with close-to-light-like separations [63, IT5]. 

The relative simplicity of the Wilson line integrals for the six-point (hexagon) case 

(2) 

allowed the remainder function FQ (ui,U2,Us) to be evaluated analytically in terms 
of Goncharov polylogarithms [HI]. The expression in ref. [HI] was then simplified 
considerably, to a few lines involving classical polylogarithms Li m , using properties of the 
symbol operation which captures all of the essential analytic behavior of a multi-variable 
function [65]. These results open the door to the possibility of finding simple analytic 
results for more general kinematic configurations. For some configurations in which the 
external momenta are restricted to two space-time dimensions, two-loop analytic results 
are already available beyond six external legs [66] • There is certainly the prospect of 
more multi-loop analytic results on the horizon, for both Wilson loops and scattering 
amplitudes in planar M = 4 super- Yang-Mills theory. 

Gravity seems to be a completely separate force from gauge theory. It has a spin 
two force carrier instead of spin one, no color degrees of freedom, and a dimensionful 
coupling constant. Nevertheless, the two theories are intimately connected. The 
AdS / CFT correspondence [2] relates the two, holographically, as a weak-strong duality. 
As mentioned earlier, strong-coupling scattering of gluons has an alternate description in 
terms of strings moving in a weakly-curved five-dimensional gravitational background. 
On the other hand, there is also a weak-weak duality of some kind between gravity and 
gauge theory, directly in four dimensions: It is possible to write perturbative scattering 
amplitudes for gravitons in terms of "double copies" of gluon amplitudes. 

The original examples of such relations are due to Kawai, Lewellen and Tye 
(KLT) [67], who found them by first deriving relations between closed and open string 
theory tree amplitudes, whose low-energy limits give tree-level relations in field-theory. 
Graviton amplitudes are represented in terms of sums of products of pairs of gluon 
amplitudes with different color orderings. The KLT relations hold not only for pure- 
graviton scattering, in terms of pure-gauge theory, but also for Af = 8 supergravity 
(or any subsector of it) in terms of amplitudes for Af = 4 sYM. More recently it has 
been recognized that other double-copy formulas exist [681 EH]- The new formulas are 
simpler in some sense than the KLT relations; they involve the direct squares of gauge- 
theory components, whereas the KLT relations employ multiple permutations of full 
gauge-theory amplitudes. As reviewed in this issue by Carrasco and Johansson |31j, 
these relations, combined with generalized unitarity, have important consequences for 



Introduction to scattering amplitudes 



11 



multi-loop amplitudes in M = 8 supergravity, because they map the problem to a much 
simpler one, the corresponding multi-loop amplitudes for (non-planar) M = 4 sYM. 

In the next two sections we discuss a few of the basic tools used to organize 
and illuminate the color and kinematic structure of amplitudes, which underlay the 
discussions in many of the other articles in this special issue. 

2. Color Preliminaries 

In this section we describe some common conventions for organizing the color structure 
of SU(iV c ) gauge theory amplitudes, which are used elsewhere in this special issue. Most 
of the articles implicitly discuss color-ordered partial amplitudes, which are particularly 
convenient in the large-iV c limit in which planar diagrams dominate. Color-ordered 
amplitudes emerge from a "trace-based" color decomposition (as also reviewed in 
e.g. refs. [6j ED]). However, in some cases, in particular for describing sub leading-color 
terms in amplitudes, and color- kinematic duality relations [31] . decompositions based 
on the SU(iV c ) structure constants are more useful. 

In general, we consider two different SU(iV c ) representations for the external states: 

• The adjoint representation, for the gluon and any of its superpartners (i.e. for 
all states in M = 4 sYM). Adjoint color indices are denoted by a, 6, c, a i; . . . e 
{l,2,...,iV c 2 -l}. 

• The fundamental (defining) representation N c , and its conjugate representation N c , 
for quarks and anti-quarks respectively. Fundamental color indices are denoted 
by *i, *2 ? - - - ^ {I; 2, . . . , N c }, and anti-fundamental N c indices by Ji,J 2 ,... G 
{l,2,...,iV c }. 

The generators of SU(iV c ) in the fundamental representation are traceless hermitian 
N c x N c matrices, denoted by (T a )\. It is conventional when discussing helicity 
amplitudes to normalize the generators by Ti(T a T b ) = 5 ab , in order to avoid a 
proliferation of \/2's. 

In QCD, the group theory (color) factors in Feynman diagrams are of two types: 
(T a )l for the gluon- quark- ant iquark vertex, and the SU(iV c ) structure constants f abc 
(or products of / afec, s) for all other vertices. (It is also convenient to normalize the 
structure constants differently, using f abc = iy/2f abc , where f abc is the standard, 
textbook normalization.) The color factors can be represented diagrammatically [7T] 
using Feynman- diagram notation, as in fig. HI Lines with arrows (quarks) carry 
fundamental representation indices; curly lines (gluons) carry adjoint indices. The 
adjoint representation is also a bi-fundamental (iV c x iV c ) representation from which 
the singlet, or trace component, has been projected out, as in the second equation for 
S ab in the figure. This relation also gives rise to the SU(iV c ) Fierz identity, 

(n<? cn<? = w-^cc, (i) 

which allows one to simplify products of traces of T a 's. 
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Figure 4. Graphical representation of building blocks for SU(-/V C ) color factors. 

2.1. Trace-based color decompositions 

Color-ordering [72J is related to the 't Hooft double-line formalism [73]. We begin by 
rewriting the color factors entirely in terms of the T a generators, using the relation 

f abc = iV2f abc = Tr([T a , T b }T c ) , (2) 

which follows from the definition of the structure constants, and is depicted graph- 
ically in fig. HI After inserting eq. fl2]) repeatedly into the color factor for a 
typical Feynman diagram, one obtains a large number of traces of the generic form 
Tr(. . . T a . . .) Tr(. . . T a . . .) ... Tr(. . .). If the amplitude has external quark legs, then 
there will also be strings of T a, s terminated by fundamental indices, of the form 
(T ai . . . T" ™)^ 1 , one for each external quark-antiquark pair. 

The number of traces can be reduced considerably by repeated use of the SU(iV c ) 
Fierz identity (fl]). In the case of tree- level amplitudes with external states in the adjoint 
representation, such as n-gluon amplitudes, fig. |5] sketches how the color factors all may 
be reduced to a single trace, Ti^T^'T^ 2 ' • • -T^t™)), for some permutation a 6 S n of 
the n gluons. This reduction leads to the trace-based color decomposition for n-gluon 
tree amplitudes, 

AT({h, K at}) = g n ~ 2 £ Tr(TMD . . . t<M»>) < ree (a(l hl ), . . . , a(n hn )) . (3) 

<r£S n /Z n 

Here A^ ee is the full amplitude, with dependence on the external gluon momenta fcj, 
i = 1, 2, . . . , n, helicities hi = ±1, and adjoint indices a^. In QCD, the gauge coupling g 
is related to the strong coupling by a s = g 2 /(4ir). The partial amplitudes (or primitive 
or color-stripped amplitudes) v4^ rcc (l fi ' 1 , . . . ,n hn ) have had all the color factors removed, 
but contain all the kinematic information. Cyclic permutations of the arguments of 
a partial amplitude, denoted by Z n , leave it invariant, because the associated trace 
is invariant under these operations. However, all (n — 1)! non-cyclic permutations, or 
orderings, of the partial amplitude appear in eq. ([3]). These permutations are denoted 
by o~ G S n /Z n = S n -x. 
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Figure 5. Reduction of color factors for n-gluon tree amplitudes to a single 
trace of T a generators. 

Looking again at the way the different trace factors arise in fig. |5j one sees that the 
partial amplitude A^ ee (l, 2, . . . , n) only receives contributions from tree-level Feynman 
diagrams that can be drawn on a plane, in which the cyclic ordering of the external 
legs, 1,2, ... ,n, matches the ordering of the arguments in A^ ee . Therefore each partial 
amplitude can only have singularities in momentum invariants formed by squaring color- 
adjacent sums of momenta, such as s^+i = (ki + k i+ i) 2 , s^ i+ i^ + 2 = (k{ + k i+ i + k i+2 ) 2 , 
etc., where all indices are defined modulo n. In this way, the color decomposition (J3j) 
disentangles the kinematic complexity of the full amplitude A^ cc - 

Similarly, tree amplitudes with two external quarks and (n — 2) gluons can be 
reduced to single strings of T a matrices, 

A'riqi, 93,..-, g n -i, q n ) = g n ~ 2 E (™ 2) • • • TM«-u) * 

(T€5 rl _2 

xA^{l q ,a{2),...,a{n-l),n q ), (4) 

where we have omitted the helicity labels, and numbers without subscripts in the 
argument of A^ ee refer to gluons. In this case there are (n — 2)! terms, corresponding to 
all possible gluon orderings between the quarks. Color decompositions for amplitudes 
containing additional external quark pairs have also been described [6]. 

We note that the partial amplitudes appearing in eq. (jlj) could also be used to 
describe the scattering of two fermions with different color quantum numbers, and (n—2) 
gluons. For example, if the fermions are gluinos in the adjoint representation, then the 
color decomposition has the form of eq. ([3]), but the partial amplitudes, with two external 
fermions, are of the type given in eq. (jlj). This is easy to see if the two fermions are 
color-adjacent; essentially all one needs to do to go between the two cases is remove one 
color line running between the two gluinos. However, the other cyclic orderings can also 
be obtained from A^ ee (l q , 2, . . . , n — 1, rig), using the Kleiss-Kuijf relations discussed 
in the next subsection. This property illustrates the universality of the color-ordered 
primitive amplitudes alluded to in the introduction. 

Many articles in this issue consider loop-level amplitudes in which all of the external 
states are gluons in the adjoint representation. The general color decomposition here is 
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similar to eq. (EJ), except that multiple color traces are now possible. Roughly speaking, 
a loop of gluons carries a fundamental and an anti-fundamental index around the loop 
with it, and an external gluon can attach to either index. The number of traces at L 
loops is equal to L + 1. 

At one loop, the full color decomposition for n gluons is [74J 



A]; loop ({k l ,h l ,a l })=g n 



N c Tr(T a -w • ■ ■ T a -<">) A n; i(cr(l' 11 ), . . . , a{n hn )) 

|n/2j+l 

+ >P Tr(T° CT « • • -T a < c -v) Tr(T aCT < c > • • •P'W) 

C=2 cr£Sn/Sn;c 

x A.Xl" 1 ),...,^^)) . (5) 

Here v4 n;c are the partial amplitudes, Z n and SV^c are the subsets of S n that leave the 
corresponding single and double trace structures invariant, and [x\ is the greatest integer 
less than or equal to x. The formula is more complicated than at tree level, but only by 
the need to keep track of more distinct trace structures and their various symmetries. 

If there are rif flavors of quarks in the loop, it is easy to see that they only contribute 
to the single-trace coefficient A n; i in eq. fl5]), and with a weight nf/N c relative to the 
adjoint gluons. (Similar color decompositions are available for one-loop amplitudes with 
an external quark pair [75J.) 

When one constructs the color-summed cross section, as is usually required for 
QCD applications, the contribution of the double-trace coefficients, A n . c for c > 1, is 
suppressed by a power of 1/N 2 with respect to that of A n -i. (It also turns out that 
the A n . c are not really independent, but can be computed as a sum over permutations 
of the A n -i [22].) Thus eq. (|5]) simplifies a lot in the large--/V<. limit. The factor of 
N c = Tr(l) simply comes from those terms in which all external gluons attach to the 
same fundamental index line, leaving the trace of the identity matrix for the untouched 
line. 

Correspondingly, if we only want the leading terms in the large- N c limit at L loops, 
we can write the compact color decomposition, 

*^7i ^ ( { ^« 5 ^ i 5 ®i } ) I leading- color 

= g n - 2 (g 2 N c ) L Tr(T a<T(1) • • • T a<T(n) ) A^ (a(l hl ), . . . , a(n hn )) , (6) 

where we have dropped the index on the leading-color partial amplitude. The 
't Hooft coupling A = g 2 N c emerges naturally in this limit. The articles in this issue 
about multi-loop amplitudes in planar M = 4 sYM are generally concerned with the 
A^ in eq. (jSJ). Because each loop integration typically brings a factor of l/(47r) 2 (or 
l/(47r) 2_<: when using dimensional regularization in D = 4 — 2e), in some cases the 
precise normalizations of the A^ will differ by this factor, as well as possibly other 
conventional factors. 
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Figure 6. The Jacobi identity ([7]). 



2.2. f -based color decompositions 

There is another type of color decomposition that is actually more useful than the 
above trace-based decomposition for addressing certain issues, and particularly for 
working beyond the large- iV c approximation. We could call this approach the /-based 
decomposition, because it goes back to the original color factors built out of the structure 
constants f abc . (Now we assume that all states are in the adjoint representation.) The 
diagrammatic representation of a general /-based color factor at tree level has the 
structure of a cubic tree graph, with vertices given by / abc 's and propagators given 
by 5 a6 's. However, these color factors are not all independent, due to the color Jacobi 
identity. This identity can be written in different ways. One way, 

Jdae Jbce Jdbe Jace Jabe Jcde 

which is also depicted in fig. El corresponds to the fact that the structure constants are 
also the SU(iV c ) generators in the adjoint representation, f abc = (F b ) ac , so that their 
commutator gives f abe contracted with F e . 

At tree level, one can use the identity (J7J repeatedly to transform all /-based color 
factors into a "multi-peripheral" form in which two selected gluons, say 1 and n, are 
always at the end of a long chain of structure constants, and the other (n — 2) gluons 
are emitted from along the ladder [76j [77] , 

A tTee ({ai}) = g n ~ 2 ^2 f aia "( 2 ) xi f Xia "(3) x * . . . j x n-3a a{n _ 1) a n 
(T£S n -2 

xAT(l,a(2),...,a(n-l),n) 
= g n-2 £ (FM 2 )...F a ^)) aia X ree (l^(2),---^(n-l),n). (8) 

Notice that this representation is identical in form to the decomposition for two quarks 
and (n — 2) gluons in eq. (BJ, except for the representation used for the SU(iV c ) generator 
matrices. Although eqs. @ and OH]) correspond to the cases of fundamental and 
adjoint representations, respectively, the same type of decomposition clearly holds for 
two external matter fields in an arbitrary SU(iV c ) representation. The fact that the 
coefficients of the /-based decomposition ([8]) are identical to the color-ordered partial 
amplitudes A^ ee appearing in the trace-based decomposition can be demonstrated by 
contracting both representations with Tr (T ai T a2 ■ ■ -T an ) and keeping only the leading 
terms at large N c [77]. 

Note that there are only (n — 2)! terms in the /-based color decomposition, in 
contrast to the (n — 1)! terms in the trace-based decomposition ([31). Therefore there are 
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identities between the A^ ee , which are group-theoretical in nature — generalizations of 
certain U(l) decoupling identities [6j [70]. These identities can be used to put any two 
legs next to each other, for example legs 1 and n: 

AT(h {«}, n, {/?}) = (-1)"' £ A T(h ^({«}{/3 T }), n) , (9) 

aeOP{a}{/3T} 

where {a} and {/?} are two sets of external gluons, whose union is {2,3, . . . ,n — 1}. 
Also, n/3 is the number of elements in the set {/?}, the set {(3 T } is {/?} with the ordering 
reversed, and OP{a}{(3 T } is the set of ordered permutations, or "mergings", of the 
two sets {a} and {/3 T } that preserve the ordering of the oti within {a} and of the 
within {P T }, while allowing for all possible relative orderings of the a, with respect 
to the These relations are known as the Kleiss-Kuijf relations [78] • They can be 
derived from eq. (jHJ) by expanding out the f abc factors using eq. (j2J), identifying the 
coefficients of the trace structures, and comparing them with the original trace-based 
decomposition [77] . 

In fact, it has been realized recently that there are not even (n — 2)! independent 
color-ordered amplitudes A^ ee , but only (n— 3)! [68], as reviewed in this issue by Carrasco 
and Johansson [31J. The additional linear relations between v4^ ree, s are not purely group- 
theoretical in nature, but also involve kinematical factors. They allow one to put any 
three external gluon legs next to each other. For example, for legs 1, 2 and n, the 
identities take the form, 

< re >(l),a(2),...,aH)= £ K^A t ™(l,2,p(3),...,p(n- l),n), (10) 

where a and p are permutations and are kinematic-dependent (but not state- 

dependent) coefficients [6H1 El]. These relations were first proved using open string 
theory [79], using a contour deformation argument [80] similar to that used in deriving 
the KLT relations |67J between gravity and gauge theory amplitudes. They have also 
been proven directly in field theory, using a recursive argument [8T] . 

Although eq. (Tl0|) expresses the consequences of certain color-kinematic identities 
in the trace basis, the identities themselves are more naturally stated in an /-based 
approach. For this purpose it is useful to not transform the /-based color factors into 
multi-peripheral form, but leave the color decomposition free, in terms of the set of all 
cubic graphs T 3 , writing 

The color factors Cj are products of structure constants for each vertex Vj, of the i th 
graph, contracted together along the propagators, Cj = Yljev, f a ^ aj 2 a ^. The factors of 
1/Pj are scalar propagators, one for each internal line in the graph, while the Ni are 
kinematical numerator factors. Consider a triplet of graphs in T^, which are identical 
except for a region from which four lines emanate. Within this region they differ 
according to the three ways of joining four legs with two cubic vertices, shown in fig. [6j 
Call the three graphs s, t and u, and let the color factors C s ^ u be normalized (signed) 
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so that they obey the Jacobi identity ([7|) as C s — C t = C u . Then the statement of 
color-kinematic duality (68] is that the associated numerator factors should be related 
by N s — N t = N u , for every such triplet of graphs. A simple and intuitive argument for 
these relations has been provided based on the heterotic string [82J. 

The existence of such a representation with local Ni (i.e. polynomials in the 
momenta) has been checked in various examples, not only at tree level, but also 
through three loops for Af = 4 sYM amplitudes [69]. As reviewed here by Carrasco and 
Johansson [31], the existence of a representation for gauge theory amplitudes satisfying 
N s — N t = N u for each triplet of Jacobi-related cubic graphs has important implications. 
Large classes of numerator factors are related to each other, and gravitational amplitudes 
are easily constructed from gauge theory ones. 

3. Kinematic Preliminaries 

Now that we have described two different ways to organize the color quantum numbers 
for scattering amplitudes, let us turn to kinematical issues, including convenient choices 
of external states and kinematical variables. 

3.1. Spinor-helicity formalism 

In QCD, the helicities of massless quarks are conserved by their chirality-preserving 
interactions with gluons. Hence it is natural to use a helicity basis for the quark 
amplitudes. In four-component notation, the spinors for the external states are taken to 
be u±(k) = 7|(1 ±7 5 )m(/c) for a quark with momentum k, and v T {k) = |(1 =F7s)f (k) for 
an anti-quark. However, in the massless limit these spinors can be chosen to be equal 
to each other, u±(k) = v T (k), and we can use two-component Weyl spinors as well. We 
will see that helicity amplitudes for external gluons can be built from the same objects. 

We want to consider amplitudes with n different momenta k iy i = 1,2, ... ,n, and 
to do so in a way that respects crossing symmetry. Thus we take all the momenta to be 
outgoing, so that momentum conservation reads yj, fcf = 0. In a realistic process with 
some incoming particles, the physical momentum for each of the incoming particles is 
simply the negative of the corresponding k^. Also, the physical helicity hi is the negative 
of the one by which we label the amplitude. (The spin S{ does not reverse under crossing, 
but hi = Si ■ ki does.) The two-component spinor for the i th state, if it is an outgoing 
fermion with helicity ±|, has several notations [51 [7U| 125]: 




(12) 
(13) 




Similarly, the conjugate spinors are 




(14) 
(15) 
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Two-component indices are raised and lowered with the two-dimensional antisymmetric 
tensors e a ^, e a/3 , etc. The spinor inner products come from contracting spinors for 
different momenta with these tensors, 



tii) = ^(xM^ih = (r\i + ) = u-ik^u+ik), (i6) 

= e^CXjU\i)p = (j + \r) = Mk-)u4k). (17) 
Spinor products are antisymmetric under exchange of labels, (j I) = — (I j), [j I] — — [Ij], 
and (j j) = [j j] = 0. Lorentz vectors can be written as bi-spinors, or 2 x 2 matrices, 
by contracting them with the Pauli matrices. A massless vector written in this way 
factorizes into the product of two massless spinors, as 



(k)aa = Ki^^aa = U + (ki)u + (ki) = (A i ) a (Aj) d , . (18) 

This factorization plays a role in the BCFW complex- momentum shift, which is 
best described in terms of the A and A variables [IB], as is done elsewhere in this 
issue [HI [251 121]. 

Amplitudes with external gluons can also be described in terms of the Aj and Aj 
variables, thanks to the spinor-helicity formalism [83]. The polarization vector s ± {k) 
for an outgoing massless vector particle with momentum k and helicity h = ±1, is 
required to be transverse to k, e ± ■ k = 0. In the spinor-helicity formalism, e ± is 
also chosen to be transverse to another massless momentum q, called the reference 
momentum (which should not be parallel to k but is otherwise arbitrary). Note that 
k and q span a two-dimensional subspace of four-dimensional space-time. Because 
|^|fc ± ) = ^|g ± ) = 0, the two orthogonal directions to this subspace are spanned by 
(g + |7 A1 |/c + ) and {q~\^(^\k~). The spinor-helicity polarization vectors live in this subspace; 
indeed, they are proportional to these two vectors, and are normalized by the condition 
that e ± ■ (e ± )* = -1: 

^>-SS?- (19) 

As a bi-spinor, e ± (k i , qi) is given by 

[e + {k t ^)U = V2 {Xq j a{Xt) \ [e-{h, qi )U = -V2 {K y X ? )a ■ (20) 

If each qi is chosen to be another momentum in the process, say kj, then it is clear from 
the general form of the Feynman rules (using also eq. ( fT8l) ) that the full amplitude can 
be built entirely out of the spinor products (j I) and [j I], for 1 < j, I < n. 

What are the basic properties of the spinor products which make them so convenient 
variables for helicity amplitudes? First of all, for real momenta, the two types of spinor 
products are complex conjugates of each other, {j I) = ±[lj]*. Also, for either real or 
complex momenta they "square" to give the ordinary dot products of momenta, 

= Tr[i(l- l5 )Mj] = 2k r h = Sjl . (21) 
(Equation (T2"T|) is written in the convention usually used in the QCD literature, but 
the reader should be aware that another convention exists, in which the sign of [j 1} is 
reversed so that {I j) [j 1} = —Sji.) 
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Figure 7. (a) A typical collinear limit in a theory with massless scalars. There 
is no suppression of the singularity from angular-momentum conservation, (b) 
A typical collinear limit in a massless gauge theory, for which there is always a 
mismatch of at least one unit of angular momentum, lessening the singularity. 



For real momenta, the conjugation property and eq. (|2"T|) imply that 

(jl) = V^Tie^ 1 , [jl] = v^e-**, (22) 

where <pji is a phase. This phase shifts as kj and ki are rotated in azimuthal angle 
around the axis corresponding to their sum, kp = kj + k\. One way of seeing the utility 
of the spinor products for helicity amplitudes is to examine the limits in which two 
particles become collinear, kp = Sji — > 0. Figure C^a) shows that the typical behavior 
of amplitudes in a massless scalar field theory is ~ 1/fcp = 1/sji, coming just from the 
scalar propagator. However, in massless gauge theory, there are numerator factors in 
the Feynman diagrams which lessen this divergence. More physically, the factorization 
in the collinear limit should be onto a physical gluon state with helicity ±1. However, 
the sum of the two external helicities must be ±2 or 0. Therefore there is at least a ±1 
mismatch of the spin angular momentum along the kp axis, as illustrated in fig. Efb). 
The spin angular-momentum mismatch requires some orbital angular momentum, which 
in turn causes a suppression of the magnitude of the amplitude, from 1/sji to 1/ \fsji- 
It also dictates a phase shift under azimuthal rotation of kj and ki about the kp axis. 
Both of these properties are captured by the spinor products, making them ideal for 
describing helicity amplitudes. 

The simplest nonvanishing tree-level gluon amplitudes are the Parke- Taylor or MHV 
amplitudes [I], 

^ r "^' pff'-CD - (23) 

Here exactly two gluons, I and m, have negative helicity, and the remaining (n — 2) 
gluons have positive helicity. The amplitudes with zero or one negative gluon helicity 
vanish by supersymmetry Ward identities |S] . The only denominator factors in these 
amplitudes are the spinor products for color- adjacent pairs of momenta, (i, i + 1). They 
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capture such universal collinear limits [SJ 122] as, 

AT e (1+ ...,/",..., m" . . . , (n - 1) + , n+) | n _ X || n (24) 
1 A t „ rec 1 (l + ,...,r,...,m-,...,P + ), (25) 



^(1 - z)(n - 1, n) 

where legs (n — 1) and n are becoming parallel, with kp = k n -\ + k n , k n -\ ~ zfcp and 
fc n « (1 - z)k P . 

3.2. Three-point amplitudes and complex kinematics 

Another way to see why the spinor products are so useful is to consider the special case 
n = 3 for the MHV amplitudes, 

(12) 4 

< ec (r, 2~ 3 + ) = % , r-V-Vi ; • (26) 
3 v ' ' ; (12) (2 3) (31) y J 

One might think that the scattering of three massless particles is singular, because 
all of the Sji vanish since they are equal to the squared momentum of the remaining 
particle. Indeed, for real momenta, eq. ( I2"2"j) implies that all the spinor products vanish 
too, and there is no way to build a nonvanishing and nonsingular amplitude. However, 
for complex momenta, the complex-conjugation relation (j I) = ±[/j]* no longer holds, 
and it is possible for half the spinor products to be nonvanishing, while the other half 
vanish. 

Specifically, we can choose the three negative-helicity two-component spinors to be 
proportional [TT] . 

Aj oc oc Af. (27) 

Then according to eq. ( fl71) we have [1 2] = [2 3] = [1 3] = 0. However, the other three 
spinor products, (12), (2 3) and (13), are allowed to be nonzero. (This is consistent 
with the vanishing of eq. ( 12T|) . which still holds for complex momenta.) For this 
choice of complex kinematics, the MHV three-point amplitude (126]) is nonvanishing 
and nonsingular, while the conjugate "MHV" amplitude, 

^< 1+ ' 2+ ' 3 ^- ! p&i]' (28) 

vanishes. In contrast, for the conjugate kinematics with 

A" oc oc A3, (29) 

eq. (1261) vanishes while eq. (1281) is nonvanishing and nonsingular. 

The three-point amplitudes ( 1261) and ( |28|) . as well as related amplitudes with two 
matter particles and one gluon, are very important in massless gauge theory. The 
BCFW relations build all tree-level scattering amplitudes recursively from them, using 
factorization onto multi-particle poles, as well as onto complex-momentum versions of 
the collinear poles shown in fig. [3 Of course Feynman diagrams can also be used to 
build amplitudes from three- and four-point vertices. However, Feynman vertices are 
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typically evaluated with off-shell lines emanating from them, which makes them gauge- 
dependent. In contrast, eqs. fl26l) and (128]) are on shell, albeit for complex momenta, 
and therefore they are fully gauge invariant. In the recursive construction of n-gluon 
amplitudes, the four-gluon vertex is never needed, essentially because it is related to the 
three-vertex by gauge transformations. 

Another reason why the three-point amplitudes are important is for building up 
loop amplitudes through generalized unitarity, as reviewed in this issue by Ita [25] and 
Britto [21] at one loop, and by Bern and Huang |30j and Carrasco and Johansson [31] 
at the multi-loop level. Typically, generalized unitarity cuts place many propagators on 
shell, so they very often contain at least one three-point tree amplitude. Finally, the 
relations between gravity and gauge theory tree amplitudes are the simplest of all for 
the three-point case, where they can be written as an exact square, 

M 3 tree (l,2,3) = [Af c (l,2,3)] 2 , (30) 

ignoring overall coupling-constant factors. 



3.3. Kinematic variables for planar theories 

For planar (large N c ) gauge theory amplitudes — which includes all tree-level amplitudes 
- new sets of kinematic variables have proved very useful for identifying new structures 
and symmetries. 

Because the external legs have a definite cyclic order in the planar case, they carve 
the plane into sectors. To each such sector, one can assign a vertex for a dual graph, 
as shown in fig. [8]^a). Lines connecting vertices of the dual graph cross lines of the 
original graph. The latter lines carry momenta. Differences between Xj's are computed 
according to the net momentum flowing through the lines of the original graph. Thus 
in fig. [H^a), adjacent sectors are separated by Xj + \ — Xj = kj, or more generally, 

x^ j = x^ Xj = kj ~\~ kj^_^ ~\~ ' ' ' ~\~ k^ ^ . (^-0 

The Xi are often called dual coordinates (or sometimes region or sector variables), but 
they are essentially momenta, not coordinates. In a planar loop graph we can assign 
additional Xi variables to the faces inside each loop, corresponding to the independent 
loop momenta. Figure E^b) shows a two-loop example. Here the thick lines of the dual 
graph that are shown each cross (and therefore correspond to) one propagator of the 
planar double box integral. 

More precisely, differences of the Xi are momenta. The Xj themselves are not 
constrained by momentum conservation. They are constrained by on-shell conditions, 

^ 3+1,3 ~ ft i 



a^+i j = k? = 0. However, unlike momenta, it is possible to invert the x^'s, according to 



Remarkably, this inversion is a symmetry of both integrands and amplitudes in planar 
M = 4 sYM. The integrands are functions of the combinations x\^ which are manifestly 



Introduction to scattering amplitudes 22 




(a) (b) 

Figure 8. (a) Dual coordinates Xj for a color-ordered tree amplitude with 
external momenta k{. The dual variables live in the sectors, or regions, 
demarcated by the ki . Thin (black) lines denote ordinary momenta, while thick 
(red) lines denote differences Xij between dual coordinates. These differences 
also equal the momenta carried by the lines they cross, (b) Dual coordinates 
for a particular integral entering the planar two-loop four-gluon amplitude. 
Now additional dual coordinates, X5 and xq, live within the loops. In this case 
we only show the seven dual coordinate separations associated with the seven 
propagators for this integral (2:1,5, £2,5, £3,5 j £1,6, £3,6, £4,6 and 2:5,6) • 



invariant under Lorentz transformations and translations of the x^s. (Translation 
invariance is guaranteed simply because the amplitude depends on momenta, which 
are differences of Xj's.) Therefore the inversion ( )32|) combines with Poincare invariance 
to generate the conformal group SO (4,2). This group acts on the dual variables Xi, 
not on space-time coordinates; hence the symmetry is referred to as dual conformal 
invariance. 

At the loop level, dual conformal invariance can be spoiled by the infrared regulator 
needed for on-shell amplitudes. A loop- integration measure in four dimensions, d A Xi, 
transforms under an inversion ( )32l) as, d A Xi — > d 4 Xi/ (xf) 4 . In D = 4 this factor precisely 
cancels terms from inverting the associated propagators in graphs such as the one in 
fig. [1(b). However, in dimensional regularization, we have D = 4 — 2e, and the invariance 
is lost due to extra factors in the transformation of the loop integration measure. As 
discussed in detail in this issue by Henn [49] , for planar Af = 4 sYM it is possible to 
retain dual conformal invariance at the loop level by using instead a Higgs regulator [48J , 
which generates particle masses that can also be thought of as extra components of the 
Xi . Although amplitude computations with the Higgs regulator have not been pushed 
quite as far yet as in dimensional regularization, there are certain advantages to this 
approach. For example, in dimensional regularization, terms in lower-loop amplitudes 
that vanish as e — > often have to be computed, because they can multiply pole 
terms, ~ l/e fe , in other amplitudes (or in the squares of amplitudes needed to construct 
differential cross sections). The article by Schabinger in this issue discusses such 0(e) 
contributions [M]. In contrast, with the Higgs regulator, the singularities have the form 
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of powers of logarithms, ~ ln fc m 2 , and power-suppressed terms of the form m 2 / 's^+i 
can always be dropped, because they vanish faster than any power of a logarithm. 

How constraining is dual conformal symmetry? In general it can only determine 
amplitudes up to arbitrary functions of the cross-ratios, 

Ui i kl = ™2 ™2 ' ( 33 ) 
x ik x jl 

because these variables are invariant under the inversion ( 1321) . However, the on-shell 
constraints xj i+1 = kf = imply that there are no non-trivial cross ratios for four- and 
five-particle scattering. At the six-point level there are three non-trivial cross ratios, 
utilizing the allowed xf i+2 and x 2 i+3 kinematic invariants: 

2 2 2 2 2 2 

X 13 X 46 S 12S45 X 24 X 51 s 23 s 56 x 35 x 62 S 34-S61 , OA \ 

Ui = 2 J 2 ° = , U 2 = 2 ° = , U 3 = f I = . (34) 

X 14 X 36 S123S345 x 25 x 41 S234S123 x 36 x 52 S345S234 

(2) 

As mentioned earlier, the remainder function R 6 characterizing the two-loop MHV 
six-point amplitude in planar M = 4 sYM is a function of these three variables, and this 
function is now known analytically in terms of polylogarithms [6U 165] . 

In planar M = 4 sYM, dual conformal invariance combines with dual M = 4 
supersymmetry to generate dual super conformal symmetry [851 ESI HZ]- It is natural 
to consider super-amplitudes, which package sets of amplitudes together by using an 
Af = 4 on-shell superfield [5], 

Hv) = 9 + + V A 9A + ^V A V B <Pab + ^V A V B V C ^abcd9 D + ^V A V B V C V D ^abcd 9~ ■ (35) 

Here g are the ±1 helicity gluons, and g the four flavors of ± 2 helicity gluinos, 
and 4>ab the six real helicity scalar states. Integrations over the Grassmann variables 
r] A , A = 1,2,3,4, can be used to pick off the desired component amplitudes from the 
super-amplitudes, defined by 

AJjh) = A($i(r]i),...,MVn))- (36) 

Just as the Xi bosonic variables automatically satisfy momentum conservation, k = 
Ya=i ki = 0, one can construct their superpartners 9f a , whose differences satisfy 

e? a -e?£ = xtvt, (37) 

for i = 1,2, ... ,n. Super-momentum conservation constitutes eight constraints: 

n 

g Aa = S^ X?rj A = Q (3g) 
i=l 

Like ordinary momentum conservation, it is satisfied by virtue of the periodicity of 

qAa 
'i+n- 



the dual coordinates, Xi = Xi +n , 9f a = 0{± n - The MHV super-amplitude generalizing 



eq. (E3D is 

,t ree ,MHV = , 5\k)5*(q) 5 \ Xl - Xn+l)5 ^ ei -e n+l ) 
n (12)(23)---(nl) (12)(23)...(nl) ' { ' 

It transforms covariantly under dual superconformal transformations [851 El] which 
extend SO(4,2) to the superalgebra PSU(2,2|4). Furthermore, the solutions to the 
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BCFW super-recursion relation for amplitudes with more negative helicities (NMHV, 
NNMHV, etc.) are given by the product of this super-amplitude with collections of dual 
superconformal invariants, generically denoted by R n [85J, EH1 El] • 

Superconformal and dual superconformal symmetry are separate symmetries, which 
partly overlap. They close into a very large symmetry group called the Yangian [88], as 
discussed in this issue by Bargheer, Beisert and Loebbert [89] and by Drummond [21J. 
The generic generator of the Yangian is not local in the sense of being a sum of single- 
particle terms (like k and q); rather it is a sum of multi-particle operators. (For the 
dual superconformal generators, the nonlocality is mild; only the sum of two-particle 
operators appears.) There are certain "anomalies" in the action of the Yangian on 
amplitudes. At tree level, the anomalies only act at the boundary of the n-particle 
phase space, where two particles become collinear, as in fig. [7] [90]. At loop-level they 
act in the full phase-space, due to singularities in the loop integration [91]. However, it 
is possible to deform either the symmetry algebra or the representation in a suitable way 
so as to preserve the invariance [901 EH [89] . It will be very interesting to see whether the 
constraints from the Yangian are sufficient to determine multi-loop amplitudes directly, 
without having to pass through the computation of loop integrals. 

Another very recent development, in the same vein of using symmetries or 
dynamical principles to bypass standard loop (or Wilson line) integrals, concerns 
polygonal Wilson loops. One considers an operator-product-like expansion for the 
Wilson loop, associated with a limit in which multiple momenta become collinear [92J. 
This expansion can be carried out at both weak and strong coupling. At weak coupling 
it has been used to compute the discontinuity of the integrated two-loop amplitude, 
directly in terms of one- loop quantities [93] . 

4. Outlook 

There has been a remarkable array of breakthroughs in our understanding of scattering 
amplitudes over the past few years, as described in the review articles in this special 
issue of Journal of Physics A. This article has represented an overview of many of 
these developments, as well as providing some basic material as an introduction to the 
remaining contributions. However, it certainly did not do justice to all of the directions 
being pursued currently. 

One of the other remarkable features of this field has been the interplay and 
cross-fertilization between formal developments and phenomenology. The Parke- Taylor 
tree amplitudes were found in the course of trying to understand patterns arising in 
the structure of five- and six-gluon tree-level amplitudes. A similar analysis of the 
structure of the five-gluon amplitude at one loop in QCD led to the development 
of the unitarity method, which has been indispensable for computing high-loop-order 
amplitudes in maximally supersymmetric gauge theory and gravity, but which has also 
produced new NLO QCD results for colliders [25J. A "radiation zero" found over thirty 
years ago in the electroweak process du — > [91] was studied for a general gauge 
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theory soon thereafter [55]; the relations found there were recognized much later as 
the four-point versions of a more general color-kinematics duality [68J. Wilson loops, 
used to characterize the effects of soft gluons in QCD, also turned out to be in exact 
correspondence to full MHV amplitudes in planar M = 4 super- Yang-Mills theory. 

This special issue should be of interest, not only to experienced practitioners in the 
field, but also to newcomers who want to get started. As the history of the field shows, 
new developments can and will come from totally unanticipated angles. Many of the 
new ideas in the future may well come to physicists whose interest in the remarkable 
aspects of scattering amplitudes was sparked, at least in part, by reading the articles 
assembled in this special issue. 
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